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Abstract 

Starting from a precise two-nucleon potential, we use the method of unitary transformations 
to construct an effective potential that involves only momenta less than a given maximal 
value. We describe this method for an S-wave potential of the Malfliet-Tjon type. It is 
demonstrated that the bound and scattering state spectrum calculated within the effective 
theory agrees exactly with the one based on the original potential. This might open an avenue 
for the construction of effective chiral few-nucleon forces and for a systematic treatment of 
relativistic effects in few-body systems. 
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1. Chiral perturbation theory for two and more nucleons became a subject of a great research 
interest in the past few years, see e.g. the pioneering work in jy, g|. One hopes to be able to 
clarify the structure of nuclear forces in this way. However, only the low-momentum matrix 
elements of nuclear forces may be systematically treated in this approach since it is based on a 
consistent power counting of small momenta and pion masses compared to the typical hadronic 
scale of Ah a d — 1 GeV (for some recent work along these lines see e.g. refs. [||, Q). A natural 
problem arises due to the appearance of shallow nuclear bound states indicating a breakdown of 
perturbation theory. Furthermore, solving the Lippmann-Schwinger equation for constructing 
the deuteron necessarily involves momenta \p\ > Ahad- For such momenta, the chiral effective 
potential constructed according to the conventional power counting rules is no longer applicable. 
This is witnessed by the fact that in the calculations for the two-nucleon system in Q, an 
additional ad hoc cut-off to tame the high-momentum components had to be introduced. To 
be more precise, this cut-off function is not commensurate with the underlying chiral power 
counting since it introduces an infinite string of local operators with increasing dimension. One 
might therefore question the validity or usefulness of such an approach alltogether. For recent 
discussions of this subject see || and a different power counting scheme has been presented 
in H . A similar problem arises in standard few- and many-body calculations based on realistic 
nucleon-nucleon (NN) forces. The potentials, if derived from meson-exchange diagrams, are 
generally based on a non-relativistic expansion in powers of momenta over the nucleon mass 
and are then used in various types of bound state equations. These usually involve integrations 
over a much larger range of momenta as used in the construction of the NN potentials. The 
same is of course also true for the various phenomenological NN forces, which are chosen more 
or less ad hoc (with the exception of the pion tail). This affects in a non-trivial way the 
calculation of observables, such as masses, levels or electromagnetic response functions. For a 
recent discussion, see e.g. ref.|Q]. In this general context the question of the existence and the 
properties of a low-momentum effective theory for nucleons are thus of great importance. We 
show here that it is indeed possible to construct an effective two-nucleon potential from a given 
realistic potential which involves only low momenta, i.e. momenta below a chosen momentum 
cut-off, but which gives exactly the same results for bound and scattering states. The cut-off 
scale introduced in our approach should be considered a physical quantity since it defines the 
Hilbert space in which the theory operates. This is different from the cut-off in a form factor 
or vertex function. It is important to stress that our approach differs from the treatment of 
the Schrodinger equation in an effective field theory framework proposed by Lepage ||. In his 
approach, an effective field theory for nucleons only is constructed for very low momenta and 
eventually pions are added. While that is certainly a valid framework, we intend to stay closer 
to the already existing nuclear physics knowledge in that our approach will eventually allow to 
match the low-momentum theory in a well-defined way to the highly successful meson-exchange 
pictures of the nucleon-nucleon force. The results presented here should therefore be considered 
as a first step in a bigger program. Finally, we remark that while it seems to be known that such 
an exact momentum space projection can be done, to our knowledge this program has never 
been carried out before. 

2. To be specific, we consider a momentum-space Hamiltonian for the two-nucleon system of 
the form 

H(p,p') = Ho(p) 6(p-p') + V(p,p') , (1) 

where 7i$ stands for the kinetic energy and the explicit form of the NN potential will be specified 
later. For illustrative purposes we stick here to a simple S-wave potential. Note, however, that 
the inclusion of spin and isospin dependent potentials can be handled along the lines outlined 
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here. Our aim is to decouple the low and high momentum components of this two-nucleon 
potential using the method of unitary transformation |], [TTJ. For achieving that, we introduce 
the projection operators 

V = J d 3 p\p)(p \ , \p \ < A , 

X = [ d 3 p\p)(p\ , \p\ > A , (2) 



where A is a momentum cut-off whose value will be specified later and r/ (A) is a projection 
operator onto low (high) momentum states with rf = 77, A 2 = A, rjX = At/ = and A + 77 = 1. 
To be precise, the separation into low and high momentum components is to be understood in 
a limiting sense, we always consider lim^oC^- — e). In this basis, the Schrodinger equation takes 
the form 

' rfHr] rfHX \ I 77* 
XHr] XHX I I A$ 




(3) 



We now perform a unitary transformation of the type 

n^n' = u ] nu , (4) 

so that rfH'X = XTi'rj = 0. The corresponding unitary operator U is parametrized in terms of 
an operator A, following Okubo 

/ (1 + AU)-V2 _At(l + A4t)-V2 \ 

\A{l + A^A)- 1 / 2 (l + Arf)- 1 / 2 j U 

and A satisfies the condition A = XAij. The requirement of decoupling the two spaces leads to 
the following nonlinear integral equation 

A (H - [A, H] - AHA) 77 = (6) 

for the operator A. In the context of the nuclear many-body theory one often introduces a 
mean field single particle basis, which defines a complete set of N-particle states. A low-energy 
subgroup of states form a model space and one is interested in effective interactions acting in 
that model space such that the same low energy spectrum results as for the underlying N- 
body Hamiltonian. A way to arrive at that effective interaction is to decouple by a suitable 
transformation the two spaces (model space and the rest space), which leads to a decoupling 
equation of exactly the form Eq.(^). In that context it is often reformulated into a linear form 
on a two-body level using the exactly known interacting two-body states (some references are 
e.g. [11] [12]). This is indeed a feasible way to proceed also in our context, as will be shown in a 
forthcoming article. Here, however, we solve directly the nonlinear equation (||). If we denote 
by q (p) a momentum from the 77 (A)-space, Eq.(|6|) takes the formp^ 

V(p, q) - J d 3 q' A(p, q')V{q', q) + J d 3 p' V(p, p')A(p' , q) 

- J d 3 q' d 3 p' A(p, q')V(q', p')A(p', q) = (E ? - E p ) A(p, q) . (7) 



#4 Our notation is such that Q(p,q) stands for the corresponding matrix element (p\Q\q) for any operator Q. 
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The quantities E^^ are the kinetic energies related to the corresponding three-momenta. This 
equation can only be solved numerically. This is most easily done by iteration starting with 

After four iterations, we then perform an average over the values of the operator A with different 
weight factors. This allows to speed up the convergence considerably (details on this procedure 
will be published elsewhere). We also provide a regularization scheme for the singularities of 
the operator A, which arise by solving this equation, i.e. at the cut-off momentum (as becomes 
obvious from Eq.(||)).[^] To be specific, we redefine the original potential V(p,p') by multiplying 
it with some smooth functions f(p ) and f(p') which are zero in some neighborhood of the point 
\p | = A and one elsewhere. The precise form of this procedure is of no interest for the following 
and will thus not be discussed in detail here. We only add that the regularization is chosen 
mild enough that it has no effect on the observables, as will be illustrated later on for a specific 
example. 

3. We now restrict ourselves to the NN S-waves. To be specific, consider a momentum-space 
Malfliet-Tjon [|l3| potential with an attractive and a repulsive part 



with t = (<fi — q-i) ■ We choose the parameters as given in [14], Vr = 7.29, Va = 3.18, 
Hr = 614 MeV and \xa = 306 MeV. Prom here on, we only consider the S-wave part of this 
potential which can be obtained analytically. Although this potential is quite simple, it captures 
essential features of the NN interaction, in particular, it supports exactly one bound state at 
E = —2.23 MeV. Since we are interested in an effective theory with small momenta only, we 
set the cut-off A = 400 MeV (or smaller). In Fig. 1, we compare the original potential with the 
effective one. The latter is defined via 

V cS = H'-H . (10) 

In the range of the small momenta, the potentials are very similar. However, one finds significant 
differences between the effective and the original potential when the cut-off, above which the 
nucleonic momenta are integrated out, is chosen very small, A < 200 MeV. This is shown in 
Fig. 2. 

4. We now consider observables. Phase shifts can be derived from the S-matrix, or equivalently, 
from a K-matrix approach. Symbolically, the relation between the S- and the K-matrix can be 
expressed as 

b ~l + iirqK' 

with q the on-shell relative momentum of the two nucleons corresponding to the laboratory 
energy being considered. We work in the framework of the latter because the K-matrix is 
purely real. This is, however, just a matter of convenience. 

The low-energy phase shifts and the bound-state energy are reproduced to a very high 
precision with the resulting effective potential acting only in the low momentum components. 
This is shown for the S-wave phase in Fig. 3 for two values of the cut-off A = 200 and 400 MeV, 
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In contrast, the effective potential is well-behaved. 
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respectively. The phase shifts from the original potential are exactly reproduced in the approach 
based on the effective potential, as long as one stays below the chosen cut-off. This is the reason 
why the solid and dashed lines in the figure fall on top of each other. The Lippmann-Schwinger 
equation in the effective approach involves by construction only momenta below the cut-off and 
thus the bound state can be calculated completely consistently. We find that the bound state 
energy is also exactly reproduced. Furthermore, the deuteron state evaluated in the effective 
theory is of course unitarily transformed. Matrix-elements of an arbitrary operator O remain 
unchanged under this unitary transformation, 

($d\0\$ d ) = (& D \0'\<!>' D ) , (12) 

with O' = U^OU, \& D ) = U^\$ D ) and $ D (p) = {p\®d} denotes the deuteron wave function 
in momentum space. For illustration, we compare in Fig. || the original and the unitarily 
transformed deuteron momentum-space wave functions for A = 200 MeV. Note that due to 
the regularization, the "spike" close to the cut-off is infact a smooth function and does not 
introduce any singular derivatives. For a larger cut-off value, say A = 400 MeV, the two curves 
fall onto of each other apart from a small interval in the vicinity of the cut-off. We thus do 
not show that case here. Note, however, that the unprojected deuteron wave function is still 
not negligible at momenta of about 800 MeV. As an illustration, we consider the expectation 
value of the modulus of the momentum operator in the S-wave deuteron, (&d(p)\P\&d(p)), 
with p = \p\. The results are listed in table p] for the two cut-off values A = 200, 400 MeV. 
Of course, for the full MT-potential we do not need the regularization. However, to illustrate 
its influence, we have also performed a calculation with a MT-potential subject to the same 
regularization as used for the effective potential (labelled "regularized" in the table). For the 
lower cut-off, the few permille deviation between the exact result and the one based on the 
regularized potential is simply due to the fact that we did not optimize the numerical solution 
of the integral equation Eq.@ to determine the operator A. If needed, one can improve these 
numbers to agree to arbitrary precision (which is not of relevance here). 





A = 200 MeV A = 400 MeV 


Exact Potential, not regularized 
Exact Potential, regularized 
Effective Potentail, regularized 


80.23 80.2308 
79.90 80.2303 
79.90 80.2304 



Table 1: Expectation value of \p\ [MeV] in the S-wave deuteron based on the exact and the 
effective MT-potential. For the exact case, we show the results with and without regularization 
at the singularity \p\ = A. 

5. In summary, we have shown how to construct an effective low energy theory for nucleons based 
on the method of unitary transformations. For a simple S-wave potential, we have shown that 
the theory projected onto the subspace of momenta below a given cut-off reproduces exactly the 
features of the original one. We hope that this study might be useful for derivation of NN-forces 
based on chiral Lagrangians in the low-momentum regime. It should also provide new insights 
into a consistent and convergent treatment of relativistic effects in few- and many-nucleon 
systems. 
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Figures 
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Figure 1: Effective two-nucieon potential (green hatched area with solid lines) in comparison with the 
original potential, Eq. (0) (blue hatched area with dashed lines) , for momenta less than 400 MeV. 
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Figure 2: Effective two-nucieon potential (green hatched area with solid lines) in comparison with the 
original potential, Eq. (^) (blue hatched area with dashed lines) , for momenta less than 200 MeV. 
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Figure 3: Phase shifts from the effective potential (solid lines) and the original potential (dashed lines) 
as a function of the kinetic energy in the lab frame. Upper (lower) panel: A = 400 (200) MeV. 
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Figure 4: Deuteron wave function p§d(p) versus the momentum p from the effective potential (solid 
line) and the original potential (dashed line) for A = 200 MeV. 
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